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Abstract 

We consider a quantum mechanical black hole model introduced in Phys.Rev., 
D57, 1118 (1998) that consists of the selfgravitating dust shell. The 
Schroedinger equation for this model is a finite difference equation with the 
shift of the argument along the imaginary axis. Solving this equation in quasi¬ 
classical limit in complex domain leads to quantization conditions that define 
discrete quasiclassical mass spectrum. One of the quantization conditions is 
Bohr-Sommerfeld condition for the bound motion of the shell. The other 
comes from the requirement that the wave function is unambiguously defined 
on the Riemannian surface on which the coefficients of Schroedinger equation 
are regular. The second quantization condition remains valid for the unbound 
motion of the shell as well, and in the case of a collapsing null-dust shell leads 

to m Vk spectrum. 
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INTRODUCTION. 


The black hole physics gives us an example of the strong gravitational fields. The exis¬ 
tence of the event (apparent) horizons causes the Hawking’s evaporation of the black holes . 
The fate of the evaporating black holes becomes a subject of interest. The quantum theory 
may throw some light on many problems of the classical black hole physics. In the absence 
of the theory of quantum gravity we have to construct a new theory every time we want to 
quantize some classical gravitating object. 

The simplest black holes models which nevertheless capture some important properties 
of the system consisting of the gravitational field interacting with matter are of special 
importance in the theory of black holes. Many authors have considered the models where the 
interaction of the gravitational field with matter in the form of a thin spherically symmetric 
shell was taken in account in the selfconsistent way ill. 

In |I| the classical geometrodynamics for the system of a thin dust matter shell in its 
own gravitational field was constructed. The canonical transformation was found which 
reduce the system of constraints to a rather simple system of functions on the phase space. 
The constraints are easily realized on the quantum level. Quantization of such a model in 
the coordinate representation leads to the Schroedinger equation in finite differences. The 
shift in the argument is along imaginary axis which has very important consequences. One 
of them is that the wave functions which are the solutions to such an equation should be 
analytical functions on the appropriate Riemannian surface. 

In the ordinary quantum mechanics we are dealing with the second order differential 
equations. We demand that the solution should be at least two times differentiable. To find 
eigenfunctions and spectrum we need to specify a class of functions, usually by imposing 
appropriate boundary conditions. In the case of the finite differences operator with the 
shift along the imaginary axis we must specify a class of functions by demanding analyticity 
(except in the branching points). In [jl[ it was the analyticity requirement that enabled us 
to find the mass spectrum in the limit of large black holes. The spectrum depends on two 


11 


2 





quantum numbers. This is explained by the fact that we succeeded in accounting for the 
motion of the shell in all regions of Kruskal space-time. The conhguration space variable, 
the radius R of the shell takes its values on the cross (see Fig. ^ with intersection at the 
horizon point R = Rg {Rg is the gravitational radius of the shell). R changes from 0 to Rg 
on two of the hands of the cross in T± regions of the space-time and from Rg to cx) on the 
other two hands in R± regions of Kruskal space-time. So the principal difference from the 
usual one dimensional quantum mechanics, where the particle moves along a real line, R^ 
was that the conhguration space was not even two copies of a real line but some nontrivial 
topological space (a cross) which is not even a manifold. The nontrivial topology of the 
conhguration space results in an appearance of the second quantum number. 

In this paper we consider the quasiclassical solutions of the hnite-diherence Schroedinger 
equation of the quantum black hole model considered in [|I|. The Riemannian surface on 
which the hnite diherence Schroedinger equation is dehned turns out to be a sphere with 
two punctured points which are two spatial inhnities at and regions of Kruskal 
space-time (Fig. |^). There is a nontrivial cycle on this Riemannian surface. The usual 
Bohr-Sommerfeld condition, that corresponds to the bound motion of the shell, gives one 
of the quantum numbers in the spectrum. Apart from this condition there appears to 
exist another nontrivial quantization condition which stems from the requirement for the 
wave function to be a regular function on the relevant Riemannian surface. It therefore 
should have trivial monodromy along the nontrivial cycle on the Riemannian surface. This 
requirement leads to another Bohr-Sommerfeld-like condition on the mass spectrum. So the 
quasiclassical mass spectrum depends on two quantum numbers. It becomes clear in our 
considerations that the appearance of the second quantum number is due to the nontrivial 
topology of the classical conhguration space of the model. 

We will remind the black hole model considered in []^ in sections | and and write 
the hnite diherence Schroedinger equation which we will analyze (it is slightly diherent 
from the one considered in because of another factor ordering). Then the method of 
construction of quasiclassical solutions of diherential equations in complex domain will be 
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explained in section El We show how Bohr-Sommerfeld quantization condition arises from 
gluing the global quasiclassical solution from solutions dehned in different regions of the 
complex plane. In sections 0 and [V| we apply this method for the Schroedinger equation of 
our black hole model and obtain the quasiclassical mass spectrum which posses the special 
properties mentioned above. It is different from the large black hole spectrum obtained in 
but coincides with it in a certain limit. In section ^ we show that one of the quantization 
conditions remains valid for the unbound motion of the shell and therefore dehnes the 
spectrum of the states of the collapsing shell. In the case of a null-dust shell when the 
continuous parameter M - bare mass of the shell is absent, this quantization condition leads 
to discrete mass spectrum for unbound motion. This spectrum turns out to be m 
k & Z - the spectrum found by Bekenstein and Mukhanov . 

I. HAMILTONIAN DYNAMICS OF THE SELFGRAVITATING THIN SHELL. 

We consider the following model of the black hole |l| . This is a self-gravitating spherically 
symmetric dust thin shell, endowed with a bare mass M. The whole space-time is divided 
into three different regions: the inner part ( I ), the outer part ( II ) both containing no 
matter helds and separated by a thin layer ( HI ), that contains the dust matter of the shell 
( Fig. g where the Carter-Penrose diagram of the space-time is presented.). The general 
metric of a spherically symmetric space-time is: 

ds^ = —N'^dt^ + L^{dr + N^dtY + R^{d6^ -|- sin^ 6d(j)^) (1) 

where (t, r, 9, 0) are space-time coordinates, A, N'', L, R are some functions of t and r only. 
The trajectory of the thin shell is some 3-dimensional surface S in the space-time given by 
some function f{t): = {(t, r, 0,0) : r = r{t)}. In the region I r < f — e, in the region 

II r > f -f e, and the region III is a thin layer f — e<r<f + e. We require that metric 
coefficients N, N^, L and R are continuous functions, but jump discontinuities could appear 
in their derivatives at the points of S when the limit e —0 is taken. The action functional 
for the system of the spherically symmetric gravitational held and the thin shell is 
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S = Sgr + Sshell = J 

I+II+IU 

+ {surf ace terms) — M j dr (2) 

s 

It consists of the standard Einstein-Hilbert action for the gravitational held and the matter 
part of the action that describes the thin shell of dnst. 

The complete set of degrees of freedom of onr system consists of the set of 
iV(r, f), N^{r, t), L{r, t), R{r, t) which describe the gravitational held and r{t) which describes 
the motion of the shell. The metric (|^) has the standard ADM form for 3+1 decomposi¬ 
tion of a space-time with lapse fnnction N , shift vector = {N^, 0, 0) and space metric 
hik = diag{L‘^, R^ sin^ 6). 

Snbstitntion of the expression (^) for the metric into the action (|^) gives in Hamiltonian 
formalism 

S = J (PlL + PrR -NH - N'^Hr) drdt + J itir- 

l+II s 

N (^R [R'] /{GL) + - 

(-L [Pl] - tt) dt (3) 


where PR{r),PL{r) are momenta conjngate to R{r) and L{r), vr is momentnm conjngate to 
f, hats denote the variables dehned on the shell snrface S, dots denote time derivatives d/dt 
and primes denote derivatives d/dr. Sqnare brackets denote the jnmp of a fnnction on the 
shell snrface: [M] = lim^^o {•A{r + e) — A{f — e)). The lapse N{r) and N, the shift N^{r) 
and fnnctions become the Lagrange mnltipliers as usnal in ADM formalism and P[ and 
are the constraints: 


H = 


Hr = 



( 4 ) 


The system of constraints contain two snrface constraints in addition to usnal Hamiltonian 
and momentum constraints of the ADM formalism. 
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ADM constraints: 


H = 0 


Hj. = 0 


(5) 


Shell constraints: 


Hr = Tl + L [Pl] = 0 

H = L2 = 0 


( 6 ) 


Karel Knchar |0 proposed some specific canonical transformation of the variables 
{R, Pr, L, Pl) to new canonical set {R, PR,m, Pm) in which Hamiltonian and momentnm 
constraints given by are eqnivalent to the very simple set of constraints : 


(7) 


Pr = 0 

m' = 0 

The idea is to nse the Schwarzschild ansatze for the space-time metric instead of the metric 


ds'^ = —F(R,m)dT‘^ -\ -^- -dR^ + 

F{R, m) 

R?{d9‘^ + sin^ OdcfP') 


( 8 ) 


where T,R and m are some fnnctions of (r, f) and F{R,m) = 1 — 2Gm/ R. One could 
construct a canonical transformation between {R, Pr, L, Pr) and (i?, Pr, m, Pm) so that the 
system of constraints (|]) is equivalent to the system of constraints (|^) in new variables in 
the phase space. 

In the presence of the thin shell the configuration space also contains the coordinates 
P, L and f. If we introduce the coordinates 


p = n + L [Pl] 


(9) 


Pr = ± 


m 


Pin 


RR! - GLPt, 

RR! + GLPl 


on the shell surface we could see that these coordinates turn out to be conjugate to P 
and f 0]. The set {m{r), Pm{r), R{r), pR^r), R, Pj^, f,p) gives the canonical coordinates in 
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the whole phase space of the system. One could consider the whole set of constraints 
and in the phase space If = {{R{r,t), Pfi{r,t), L{r,t), PL{r,t), r{t), n{t))}. The surface 
momentum constraint = 0 takes the form 

p=0 ( 10 ) 


The shell Hamiltonian constraint is expressed through the variables as follows 


which means that 


out 


exp 


GP, 


R 


P 


M = 0 




F, 


out 


Fir,, exp 


GP. 


R 


P 


where 


exp 



‘^\lFnFout 


Fin T Font 


M^G^ 

P^ 


± 


2MG \ 

—0 





4iH^ 


- {Font + Fin) + 




( 11 ) 


( 12 ) 


(13) 


Since the shell could be found in each of the four regions R±,T± of Kruskal space-time 
the dynamical variable R which measures the radius of the shell E embedded in the space- 
time M could take its values on the cross, shown on Fig. ^ So the conhguration space of 
the dynamical system under conhguration has the singularity at the horizon. We will see 
how this singularity disappears when we turn to the quantum mechanics where we’ll have 
to consider the realization of the relevant operators on some complex Riemannian surface 
(section 0)- 

The Hamiltonian constraint (PH|) was derived under the assumption that both Fin and 
Font are positive. It is possible to derive analogous constraints in T±-regions, where F < 0. 
But, instead, one could consider the function F^F ^ complex valued function. The point 
of the horizon F = 0 becomes a branching point , and we need the rules of the bypass. We 
assume the following 
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^ 1/2 ^ 1^1 
0 = 0 in i?+-region 

0 = 7r/2 in T_-region (14) 

0 = TT in i?_-region 
0 = —7r/2 in T+-region 

The reason for such analytical continuation is that we will be able to get the single equation 
on the wave function T which covers all four patches of the complete Penrose diagram for 
the Schwarzschild space-time. Some important consequences of this fact will become evident 
in section |V[ 

The case of special interest for us will be the dynamics of the null-dust shell which 
corresponds to the case M = 0 so that the shell propagates with the speed of light [|^. In 
this case the shell constraints take a simple form 


Hr = TT + L [Pl] = 0 

R \R'] ^ r) 

and from (PTf ) we get the form of shell constraint in terms of Kuchar variables: 

= 0 


(15) 




Fin + Fc 


out 


- (GP'^ 

2\lFi„Fnut ch ' -^ 


R 


(16) 


which is equivalent to 


exp 


GP 


R 


F 

J- 7.r 


(17) 


\ R J V Font 

In the rest of the paper we will restrict ourselves with the motions of the shell when 
rriin = 0. In this case it is convenient to make a canonical transformation from {R, to 

{s,Psy. 

ip ^ ip 

TT? d2 

(18) 


5 = 


(2Gm) 


p. = RlpA 


s - ^^g2R 

where Rg is the gravitational radius of the shell. Dimensionless variable S is the surface 
area of the shell measured in the units of the horizon area of the shell of mass m. 
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II. DIRAC QUANTIZATION OF THE MODEL. 


The Dirac quantization of the black hole model under consideration looks like the follows. 
The phase space of our model consists of coordinates (i?(r), P/j(r), m(r), 

Pm{r), R, Pj^,r,pr) r G {—oo,f — e) U(f^ + e, oo). Then the wave function in coordinate 
representation depends on configuration space coordinates: 

T = T(P(r), m(r), P, f) (19) 


and all the momenta become operators of the form 

Pnir) = -i 6/6R{r)-, Pm{r) = -i 6/6m{r)-, 

Pf^ = —i dj dR; pr = —i d/ df. (20) 

Using the Kuchar constraints (|^ and the shell constraint ([T0|) in operator form we conclude 

that the wave function does not depend on P(r) and f as far as 

/ 

d^jdRir) = 0 

< m'(r)T = 0 (21) 

9T / 9f = 0 

The dependence on m(r) is reduced in regions I and II to T = 5(m — m±) where m± do no 
depend on r. m± equal to Schwarzschild masses in the inner and outer regions min and mout- 
We restrict ourselves with the case when = 0 so the dependence of the wave function 
on the phase space variables reduces to 

\l/ = T(m, P) (22) 


The only nontrivial equation is the shell hamiltonian constraint (|^. It contains the square 
root expression which is difficult to realize at quantum level. So we use the squared version 
of the shell Hamiltonian constraint, which in terms of the canonical pair {S, Pg) (p^ reads 
as 


C 


1 

27s 


Sm'^S 


7f ch 


2Gm‘^ J 


= 0 


(23) 
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The operator C contains the exponent of the of the momentum Pg. This exponent 
becomes an operator of finite displacement when Ps is the differential operator: 

mil d 


exp 


f—1 

\2Gm^ j 


T = exp -i- 


T = T(5 - C^) 


(24) 


2m^dS, 

where mpi is Plank mass and C = 

The constraint C becomes an equation in finite differences if we substitute the expression 


(P5|) into (pdl) . 

T(m, S + iC) + T(m, S - iQ = 




\h(m, S) 


(25) 


y/S 4m^S^ 

The shift in the argument of the wave function is along an imaginary axis. In the case 
of differential equation we require the solution to be differentiable sufficiently many times. 
Similarly, we have to demand the solutions of our finite differences equation (^) to be 
analytical functions. This condition is very restrictive but unavoidable. The importance of 
this requirement is shown in where it is the analyticity of the wave functions and not the 
boundary conditions that lead to the existence of the discrete mass (energy) spectrum for 
bound states. In the next section we will see how it works in the quasiclassical regime. 

The construction of quasiclassical solutions of the differential (finite difference) equations 
in complex domain requires the use of a special technic |l^ explained in the next section. 
It is a well known fact that the quasiclassical approximation is not valid in the vicinity of 
the turning points of the classical trajectories of the system. When solving the equation in 
complex domain one could use the quasiclassical ansatze for the approximate solution only 
in regions that are distant from the turning points on the complex plane 0. One need to 
glue the global solution defined on the whole complex plane from the solutions defined in 
different regions. This global solution must be an approximation to some analytical solution 
of the differential equation under consideration. 

The next section is devoted to the explanation of the method of constructing the quasi¬ 
classical solutions of a differential equation in complex domain |OT using a simple example 
of nonrelativistic Schroedinger equation for the particle moving in a potential well. Then 
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in section ^ we use this method in order to build the analytical quasiclassical solutions of 
(p5|) in complex domain which will enable us to hnd the quasiclassical mass spectrum of our 
black hole model. 


III. QUASICLASSICAL SOLUTIONS OF SCHROEDINGER EQUATION IN 

COMPLEX DOMAIN. 


In this section we explain the construction of quasiclassical approximation to a regular 
solution of Schroedinger equation on the complex plane. The method was developed in 
and IP where all necessary theorems could be found. 

Let us consider the nonrelativistic Schroedinger equation 


- SV'(2) + ?(2)'I'(2) = 0 


(26) 


on the complex plane z E C. Let the function q{z) be holomorphic in a region D E C. Then 
in accordance with Cauchy theorem p equation (^bf) has a solution regular in D. 

We look for the approximate solutions of (^) in the form of quasiclassical ansatze 


flz) = exp ^-S)(z)j ^ ftVi 


(27) 


Then substituting ( pT]) in ( ]26| ) we obtain in zero order on h the Hamilton-Jacobi equation 
on S{z)\ 


'on'' 

dz 


q{z) = 0 


(28) 


If q{z) 7^ 0 in the neighborhood Uz^ of some point Zq then (^SD has two solutions correspond¬ 
ing to the two different branches of the function q^^'^{z): 


Z 

^2(20, ^ / J-7(2)(22 


(29) 


ZQ 


and 


= a^i{z) + b^2iz) 

4’i,2(2) = («( 2 ))^‘'’'‘exp I J ■J-q{t)dt 


( 30 ) 


^0 
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The coefficients of (^61) are analytical functions in the neighborhood of zq. Hence in ac¬ 
cordance with Cauchy theorem in the neighborhood Uq of a point q^Zo) =0 exist two regular 
solutions of (p 6 |) as well. Nevertheless, if we consider the equation (p 6 | ) far enough from zq 
where the quasiclassical approximation is valid, we hnd that the quasiclassical ansatze is 
not a good approximation for the regular solution in the punctured neighborhood of the 
point zq (the branching point for the function q^^‘^{z)). This could be easily seen if one 
note that $1 transforms into after the analytical prolongation along a closed path 71 
which goes around zq (Fig. H) while the regular solution have the trivial monodromy along 
this path. The approximate expression (|30D for the wave function is valid only in a certain 
sector a < Arg( 2 ; — zq) < (3 in the neighborhood of Zq. The explanation of the effect is the 
following. 

Let z = Zq be, for example, a simple zero of the function g : q ^ q^i^z — Zq) in Uq 
{ qo = and {z — zq) = pe*'^). Then 

nU. ^o) = 1%'Hz - (31) 

This function grows when p = \{z — Zo)\ —>• cx) in the sector —vr < + <po ~conhned 

by the lines h and I 2 ( Fig. where 0o is taken to be 0o = vr), then it is pure oscillating on 
li, I 2 and Z 5 and is decreasing in sectors vr < (30 -1- 0o ~ < 3vr and —vr < (30 -1- 0o) < —3vr. 

So the solution $1 (^) grows in the sector 1 and decreases in the sectors 2 and 3 on Fig. 

The solution $2 has an opposite behavior: it grows in sectors 2 and 3 and decreases 
in sector 1. Therefore in sector 1 the solution $1 is exponentially large compared to $ 2 - 
But in quasiclassical expression for the wave function we can not retain the exponentially 
small items simultaneously with exponentially large in the approximate expression for the 
wave function and in fact we could set the coefficient b to be equal to any number if only 
a 7 ^ 0. The explanation for the fact that the quasiclassical ansatze fails to be the correct 
approximation for the regular solution is that the coefficients a, b (^) of expansion of the 
wave function on the system of solutions $ 1^2 are different in the sectors 1,2 and 3. This 
fact was observed for the hrst time by Stokes and the phenomenon is called the Stokes 
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phenomenon [^. The lines U dehned by the equation 


Im Zq) = 0 


( 32 ) 


are called the Stokes lines. They start at the branching points q{zo) = 0 and terminate at 
singular points of equation or at inhnity. Both solutions $ 1^2 (0 are oscillating along these 
lines. 

When the point zq is a zero of n-th order for the function q{z) ( q{z) ~ qo{z — Zq)'^ ) 
the Stokes lines approach the rays Arg {z — zq) = (pk = const near zq- (pk are given by the 
formula 


(pk — 


{ 2 k + l)7r - (po 
n + 2 


(33) 


where (po = Arg (go)- The formula remains valid when n < 0 and zq is a singular point of 


The global structure of the Stokes lines when q{z) has the form of potential well with 
two turning points g( 2 ;o,i) = 0 is shown on Fig. |^. The Stokes lines divide the complex plane 
into the sectors 1,2,3 and 4 (Fig. |^) in which one of solutions (PUD is exponentially small 
compared to the other. 

Let the approximate expression for the regular solution be given in one of the sectors, 
for example in 1 


'^{z) = ai^i{z)+ bi^2{z), ze 1 


(34) 


The form of approximate solution in other sectors of the complex plane could be found using 
the using the following algorithm ||TO . 


1) The whole complex plane is divided into the set of intersecting regions {Di} called the 
canonical regions. The boundary of each region consists of some Stokes lines dDi = U4 
so that each Di is a sum of the sectors 1,2,... into which the Stokes lines divide the complex 
plane. For example, for the problem of a particle in a potential well (Fig. some of 
canonical regions are Di = 1 2 , D2 = 31J3 and D 3 = 2 IJ 4 . The regions Di 
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are foliated by the level lines Lc = { Im(r2(2;)) = c = const}. Each Di is chosen so that 
c = Im(r2) changes from —oo to +cxd in Di. In each of these regions the fnndamental system 
of solutions $ 1^2 (^) 


d'(z) = + a2<h2‘(^), zeDi 


^ 1,2 — 




exp 


(35) 


is introduced. In order to dehne the fundamental system of solutions one need to choose 
a branching point Zq G dDi and a Stokes line k G zq G k. Then the phase hi of the 
quasiclassical wave function (|^) is 


Z 

f*''( 2 i 20) = y + SJo 


(36) 


ZQ 


where hlo is chosen so that Im (hi**) = 0 on /j. Assume the hrst solution of the fundamental 
system is increasing in the sector on the right (clockwise) from the reference Stokes line and 
is decreasing in the sector on the left (anti-clockwise) from the reference Stokes line. 

2) In intersections p] Dj the transition matrices Tij between the corresponding systems 
of solutions $12 cind $ 1(2 are dehned: 


S' = aii'i + bit'i = ajt‘{ + z € 


/ 

\ 


/ \ 

di 

= T 

aj 

V 





(37) 


3) Let the form of the solution be given in some region Di. Then we could hnd the form 
of the quasiclassical solution at any point A of the complex plane. We need to draw a path 
7 from a point B ^ Di to A and to cover it by the set of intersecting canonical regions 
Di,..., Dk. The form of the solution in the region which contain the point A, is 


where the matrix T: 


T(p = +6fc<h2p z^D 


A*) 


= T 


\ bi i 


(38) 
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( 39 ) 


T — Tk^k-i X Tk-i,k -2 X ...T24 

Let us consider the quasiclassical solutions of equation (p^) in the simplest case when the 
function q{z) has two zeros (the turning points of the classical motion) situated in points zo 
and Zi on real line (see Fig. We impose on the wave function the requirement tt'(-s) —0 
when z —^ ±cx) along the real line. So the solution in the region Di = 1 IJ mbox2 must be 

^(z) =ai^[^+bi^i\ zeDi 

laA 

since the second solution of the fundamental system in region Di is decreasing in the sector 
1 in accordance with the convention taken for enumeration of the solutions from the funda¬ 
mental system. We need to find the solution in the region D 3 = mbox2 U 4 which contain 
the real line on the left from the region of the classical motion 

'F(z) = 03$!^ + 2 e l>3 (41) 




0 

1 


(40) 


We choose a path 72 connecting points A E D 3 and B E Di (Fig. ||). It could be covered 
by the three canonical regions Di, D 2 and ZI 3 so that 




fo] 


= T 


v'’’/ 




(42) 


where 


T — T3^2 X T2,2 X T24 


(43) 


The matrix T 2,2 has appeared in the equation because in the canonical region D 2 there exist 
two different fundamental systems of solutions. One of them is defined by the Stokes line I 2 
and the branching point Zq while the other is dehned by the Stokes line I 2 and the branching 
point Zi (Fig. ^). The transition matrix between the two fundamental systems of solutions 
in D 2 is 
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To , = e*' 


exp 


0 exp 

/ \/-(l{z)dz\ 


0 


^0 


(44) 


where is some factor coming from the normalization of the solutions, and which is not 
important for us. The transition matrix between the fundamental system of solutions dehned 
in the region Di with the Stokes line h with the turning point zq and the fundamental system 
of solutions dehned in the region D 2 with the Stokes line I 2 and the same turning point Zq is 


T2,i = 




(45) 


The same transition matrix is between the region D 2 with the Stokes line I 2 and the turning 
point Zi and the region with the Stokes line I 3 and the turning point Zi. In accordance 


with (|43|) we hnd 


T = 


exp 


V 


{i^{zo,Zi)) 

0 


zexp (^^n{zo, zi)'^ 

i (exp (f ^ 1 )) + exp (^-^n{zo, ^i))) 




/ 


(46) 


where fI(2(o, ^ 1 ) = I J—q{z)dz. Since we impose the requirement that the quasiclassical 
20 ^ 

wave function decreases in the sector 4, then the expansion of the wave function on the 
fundamental system of solutions in region must be 


\h) 


vV 


(47) 


because the hrst solution of the fundamental system dehned by the Stokes line /s decreases 
on the left from l^. From (^) and (^6]) we see that this requirement is satished only if 

21 


exp 


\ j \/-4(z)<*z I = 


(48) 


ZQ 


This condition is equivalent to the well known Bohr-Sommerfeld quantization condition 

21 


1 

nh 


J \/-q{^)dz = n + -, n e Z 


(49) 
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So we have found that in the simplest case when the regular function q{z) from the 
Schroedinger equation (p^) has just two simple zeros q{z) = 0 situated on the real line, 
the quasiclassical solution dehned on the whole complex plane exists if the coefficients of 
the equation satisfy Bohr-Sommerfeld condition (^91). The approximate expression for the 
solution of the form (P?] ) is different in the different canonical regions Di of the complex 
plane and the approximation itself is valid only far from the turning points q{z) = 0 of the 


equation (|26|). 


IV. QUASICLASSICAL SOLUTIONS OF SCHROEDINGER EQUATION FOR 

SELFGRAVITATING DUST SHELL 

In this section we will hnd the quasiclassical solutions in complex domain for the equation 
(p5|) which is the quantum form of Hamiltonian constraint for the system that consists of a 
selfgravitating thin dust shell and its own gravitational held. 

A principal difference of the equation p5| ) from the nonrelativistic Schroedinger equation 
(p 6 |) that was considered in the previous section is that its coefficients contain the function 

= XTT = (50) 

VP P 

(where p = y/S) that is a branching analytical function on the complex plane. So we have to 
consider the equation (pSD not on the complex plane but instead on the Riemannian surface 
Sp on which the coefficients of the equation are regular functions. 

The Riemannian surface for the function (|50|) is a two-dimensional sphere consisting 
of two Riemannian spheres 5*+ and glued along the sides of the cuts made in both spheres 
along the interval (0,1). Let us consider the real section Im(p) = 0 of the Riemannian surface 
Sp- Its part p > 0 turns out to be a cross with two ends situated in the points p = -|-cxd 
of Rf (see Fig. |^). so that > 0 in p = -|-cx) in Rf ( cx) at R+ domain ) and < 0 
in p = -l-cxo in S'_ ( cx) at F_ domain ). Two other ends are in the points p = 0 and the 
argument of F^/^ is -|-7r/2 ( T_ domain) in one of the ends and —7r/2 ( T+ domain ) in the 
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other. So we could naturally identify the cross in the real section of the Riemannian surface 
Sp with the cross discussed in the section | which is the classical conhguration space of the 
dynamical system under consideration. 

An important feature of the form of Hamiltonian constraint for the black hole model 
introduced in is that by ascribing the special choice of the argument for the function 
as a complex valued function one could obtain the constraint which is valid in all the four 
regions of Kruskal space-time and T_. The price for this was that the classical 

conhguration space of the dynamical system had a singular point at the black hole horizon 
and the conhguration space turned out to be the cross (which is not even a manifold). 
Nevertheless, when we pass to the quantum mechanical problem, we have to consider the 
complex Riemannian surface Sp as the conhguration space for the system which is a manifold 
and, as it will be shown, the dynamical system has no singular point at the horizon if we 
choose the coordinate which is regular on the surface Sp in the neighborhood of the horizon. 

We will construct the quasiclassical wave function of the equation (p5|) on the whole 
surface Sp. According to the method explained in the previous section we need hrst to hnd 
the turning and singular points of the equation. Then we have to determine the structure of 
Stokes lines and canonical regions on the Riemannian surface Sp. The last step is to hnd the 
corresponding fundamental system of solutions in each region and the transition matrices 
between them. 

In the limit of large black holes the displacement parameter ( = ‘IrrS becomes small 
and we could cut the Tailor expansion of "^{S + i() on the second term |I|. Then equation 
(p5|) becomes an ordinary differential equation: 


- + 


1^2 - F - 1/2 




T(R) = 0 


(51) 


This approximate equation is valid rather far from the horizon (|p| S> 1) in both S^ and 
S- components of Sp. So we expect that the two equations (p5| ) and (|5l|) have common 
quasiclassical solutions in the vicinities of p = cx) in 5"+ and S-. We will use this assumption 
in order to glue the quasiclassical solutions of (^5[) in different canonical regions in the 
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neighborhood of inhnities. 

One important note should be made at this point. We have seen in the previous sec¬ 
tion that the quasiclassical ansatze could not be valid in the neighborhood of any turning 
point {q{z) = 0) of nonrelativistic Schroedinger equation. Nevertheless, the Cauchy theorem 
guaranteed that there exists a regular solution of the equation in a domain containing this 
point and we could build its quasiclassical approximation everywhere except for the neigh¬ 
borhoods of the turning points. Now the picture is different. We have no theorem for the 
equation in hnite differences which affirms the existence of regular solutions of the equation 
with regular coefficients. So the assumption we made is in fact the assumption that there 
exists a solution of the equation (^) regular everywhere on Sp except for the singular points 
of the hnite difference equation (^51) (which are p = oo at Sp and S', components of S'iT’ as 
we will see) and we look for its approximate expression in the form of quasiclassical ansatze 
in the regions located far from turning and singular points of (|^) on Sp- 

For the massive dust shell there exist three qualitatively different types of classical mo¬ 
tion, described in []^. 

a) when m < M < 2m the trajectory has the form shown on Fig. ^a. It starts from p = 0 
in Tp region, crosses its own horizon p = 1, expands to some maximal radius Pmax > 1 in 
Rp region (it could be observed by an observer at inhnity during this period) and then it 
collapses to the singularity in the T_ region (the situation is called ’’the black hole case” in 
terms of |l|). 

b) when M > 2m the shell is in R_ region at the moment of its maximal expansion as it 
is shown on Fig. 0b, so for the external observer at Rp inhnity it does not appear from its 
horizon during the whole evolution (’’the wormhole case”). 

c) the case M < m describes the situation of collapse (Fig. 0c) when the trajectory goes 
from p = cxo in Rp region to p = 0 in T_ region. 

We will consider in the next section the wave function of the black hole case (a). 
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V. THE QUASICLASSICAL WAVE FUNCTION IN THE CASE OF BLACK 

HOLE. 


Let us choose the wave function in the form 


vl/(^) = exp j (HS) + H{S) + ...) 

of quasiclassical ansatze. Then the phase fl{S) satishes the Hamilton-Jacobi equation 

2 / 71 f2 \ 


(62) 


i) 


AC 


2a/5 Srri^S 


= 0 . 


(53) 


Its solutions are given by the functions O = / Ps{s)ds: 

So 


Ps = \n(F-^/^ (l + ^+ ^ 


± 


M 




where Z is given by 


2VS Stu^S 2my/S' 

In terms of the coordinate p on the Riemannian surface Sp 

P = Vs 

P, = 2pP, 


(54) 


(55) 


we have 


Pp = 2pln 




1 M ^ 

1 + v-h ::—n—T i z-Z 


2 p ' Sw?p^ 2mp 
The turning points Pg = (90/ dS) = 0 are the solutions of the equation 


(56) 


Z = 


Pi = 


1 


= 0 




Am‘^{M/ m — 1) 




P2 = 


Am^{M/ m+ 1) 


, in region 


, in I/_ region 


(67) 


one of them is in Sp component of Sp (this is the maximum radius of expansion of the shell) 
and the other is in S- component and is situated in the region p < 0 which is denoted as 
U_ on Fig. In the neighborhood of these turning points Pg is close to zero and in the 
equation 
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sh (Ps) = 


ZM 


we could set sh (P<j) ~ Pg and find that in terms of regular coordinate p 


( 68 ) 


Pf ~ \lf> - t>h ! = 1,2 


(59) 


SO the turning points are just like simple turning points of nonrelativistic Schroedinger 
equation considered in the previous section. In accordance with formula (^) we find that 
three Stokes lines originate from each of these turning points as it is shown on Fig. ^ 

In the neighborhood of infinity in region the regular coordinate is 


^ ^ 

p, = -A 


(60) 


2 P 

The solutions of Hamilton-Jacobi equation in the neighborhood of infinity are 


Pt = —^ hi I ^ _ 

2 P \ 




M 2 


± 


M 




8 m? 2m 

The argument of logarithm is close to one so taking sh (Pg) ^ Pg we obtain 


(61) 


Pt Ki ±— 


/M 2 


8P V 


(62) 


So the infinity in region is a singular point of the 4-th order and from (|33|) we find that 
the two Stokes lines originate from this point as shown on Fig. |^. This means that the 
Stokes lines in the 5+ component of the Riemannian surface Sp have vertical asymptotes. 

In the neighborhood of infinity point in R_ region Hamilton-Jacobi equation takes the 
form 


<=“ = -phi b ■ 5* 


M2 

8m2 


(63) 


and 


2PtP ^ iivr; Pt ^ ^ 


(64) 
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This singular point is of the 6-th order with four Stokes lines originating from it, as it is 
shown on Fig. ^ The Stokes lines in component of Sp approach asymptotically the lines 
tilted with the angles (j) = ±7r/4 with respect to the real line. 

The regular coordinate in the point p = 1 is 


I u = Vp-l 

[P„ = Au{l + u^)P, 

In terms of this coordinate 


(65) 


^ 1/2 _ 


U 


Vl + u^ 


( 66 ) 


and 


M 2 


M 


Pu = 4m(1 -|- u^) In < — vTT~u2- -—- , ± - 

\u \ 2a/ 1 -h m 2 8m2(l m 2)3/2 2mA/l -h u'^ 


Z]\ (67) 


The coefficients of the equation ( [25|) have singularities in the point p = 1. But we see from 
(|67|) that Pu —>• 0 as M —0 and the phase 


n = y PsdS = J Pudu = —2u'^{lnu — 1/2) -|- 

{regular part) 


( 68 ) 


also remains hnite as m —> 0. So in terms of the regular coordinate on the Riemannian 
surface Sp in the neighborhood of p = 1 this point is in fact a turning point (P^ = 0) of the 
quasiclassical solution. 


If we consider the truncated equation (|^) in terms of regular coordinate on near 
the horizon we will see that the point u = 0 is the turning point for this equation as well. 
Indeed, Hamilton-Jacobi equation for (^Tj) has the form 


P/ 


16 m 2(1 + u 


2\2 


= -{ 2 VI + 




u 


vT+ 




M 2 


4m2(l -1- ^2)3/2 


- 2 


(69) 


and 


22 


















Pu ~ a/m 


( 70 ) 


in the neighborhood of m = 0. So the turning point is just simple turning point and three 
Stokes lines originate from it as shown on Fig. |^. Note that the i?_, T+ and T_ lines 
intersect at the horizon as it is shown on Fig. So one of the Stokes lines coincides with 
i?_i_ while the two others lie between and T± lines which means that they both belong to 
S- component oi Sp- 

We are looking for the quasiclassical solution of the hnite difference equation (^5|) in 
regions far from the turning points. In these regions the approximate equation (|^) is valid. 
So we expect that the two equations have the same set of approximate solutions. So we 
will use the transition functions between the fundamental systems of solutions in different 


canonical regions (see previous section) calculated for the quasiclassical solutions of (|^ 
in order to glue the quasiclassical solutions of (^) in different canonical regions. As it 
was noted at the beginning of the section there is no general theorem about the existence 
of solutions of hnite difference equations. The method we use is based on our assumption 
about the existence of regular solutions of (pSf ) in the neighborhood of the turning points and 
they could be approximated by the quasiclassical solutions in regions far from the turning 
points. 

The remaining singular point is p = 0. Again, we pass to the regular coordinate on Sp 
near this point: 


(71) 



then in terms of this coordinate 




(72) 


and the solution of Hamilton-Jacobi equation is 
F,, = Av^ In 


\/l — 


1 m2 M ^ 


2 v 8m‘^v^ 2mv 


(73) 


23 










It follows from (^) that p = 0 is a turning point as well as the horizon because > 0 
with n —>■ 0 and the phase hi (|6^) of the wave function 


n = n"^(lnn — 1/4) + {regular part) 


(74) 


remains hnite in the vicinity of p = 0. Using the same procedure as at the horizon we 
hnd the behavior of the Stokes lines at this point using the truncated equation The 

Hamilton-Jacobi equation for (^) is 


P/ 


2 v - 


M2 


J- 


(75) 


16n® " y J 

and n = 0 is the turning point of third order. In accordance with (|33D we draw hve Stokes 
lines from this point. One of the Stokes lines coincides with U_ {Arg v = vr/2)) as far as 
the argument of logarithm in (|73D becomes zero on this line and the other four Stokes lines 
are situated symmetrically with respect to 17+ U U_ line (see Fig. |^). Two of them are in 
component of Sp and the two others are in component. 

The general structure of Stokes lines on the Riemannian surface Sp is presented on Fig. 
1^. The Stokes lines divide the Riemannian surface into eight regions. In regions 2 and 4 on 
the picture the phase of quasiclassical wave function changes in finite limits a < Im(r2) < b. 
In other regions the phase changes so that 0 < | Im(n)| < cx). Given a Stokes line and 
a canonical region containing this Stokes line we could construct a fundamental system of 
solutions in this region. 


Now we could construct a quasiclassical solution of (^5]) in the whole Sp. We will look 
for the solutions obeying the two following natural requirements: 


• (A) the solution must decrease when p —>■ cx) along the P+ and P_ lines (we consider 
the wave function of bound motion when the shell could not propagate to inhnity); 

• (B) the solution must be an unambiguous function on the Riemannian surface Sp. 

( 76 ) 
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So we start from the decreasing solution in the sector 1 (Fig. The canonical region which 
containsi?+ is = 1 U 2 U 5 with the reference Stokes line li, then the decreasing 
solution is the second solution of the fundamental system: 


T = ai^[^ + P ^ Di 


yf>i/ 


\h 


(77) 


In order to hnd the form of solution on line we should continue the solution along some 
path 7 , connecting the points A on and B on R_ (see Fig. |^). To do this we need to cover 
7 by the set of overlapping canonical regions Di, D 2 = 21J4U5U6 and = 3 U 2 U 5 (see 
Fig. |7|) and hnd the transition matrix between the fundamental system of solutions $1 and 
^2 in and <Fi, $2 in D^\ 


T 


= T ^^2 X T 2^2 X T24 
= T 






yhj 


VF 


(78) 


The matrix T 2^2 arises by the same reason as in the transition matrix (^3]). In the canonical 
region D 2 two different systems of solutions are dehned by the Stokes line I 2 and different 
turning points p = Pi and p = 1. 

The solution $ 2 ® is decreasing on R_. So the requirement (^,A) takes the form 




+ ^ 34 * 2 ® 1 p ^ Dz 


( . \ 


03 


(A 


(79) 


\h) (^7 

As far as both turning points involved in the calculation are simple turning points, the 
calculations are the same as in the previous section and we hnd that in order to obtain 
the quasiclassical solution decreasing with p —>■ 00 along both and R- lines the Bohr- 
Sommerfeld quantization condition must be satished on the Stokes line I 2 (Fig. |^: 



h 
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(80) 














Let us note the important property of the constructed solution. Consider the form of 
solution on T± lines. Let us recall the classical behavior of the shell in regions T± of space- 
time. T_ is a nonstationary region of inevitable contraction and the classical solution for 
the shell is just collapsing solutions while the T+ region is the region of inevitable expansion 
and the shell trajectory starting from singularity p = 0 expands up to p = 1 and leaves 
the T+ region (Fig. p. The quasiclassical wave function for both (^) and (^T]) equations 
expresses through the two basic solutions 


T = exp ^±—f2(S')^ (j) 


(81) 


corresponding to in-going and outgoing waves respectively. So it seems at hrst sight that at 
quantum level the shell does not contract inevitably in T_ region. But the phase Q{S) = 
f PgdS is complex on this line 


= In 


ir-i/2 1 + 


+ 




± 


M 


2VS Sm'^S 2my/S 


-.Z 


— i Arg(F^/^). 


(82) 


In each of the regions 1..8 on which the Stokes lines divide the Riemannian surface Sp one 
of the solutions (|8TD is exponentially small compared to the other. The function in the 
neighborhood of the horizon p = 1 is given by 


pl/2 ^ » 

VTT^ 

and the argument of P^/^ coincides with the argument of the regular coordinate u on Rie¬ 
mannian surface Sp in the vicinity of p = 1. The argument of u is zero on the line Rp so the 
solution is an oscillating function on this line near the horizon. On the line P_ the argument 
of u is ±7r so we have decreasing and increasing components in the fundamental system of 
solutions in sector 3 of Riemannian surface Sp. The same situation is on lines T±. The 
argument of u is Si'k/2 and the solution T = exp(-(-if2/h) which represents an outgoing wave 
is decreasing on line T_ while the solution T = exp(—zf2/h) (an in-going wave) is increasing 
and therefore is exponentially large with respect to the outgoing wave on the line T_. On 
the line T+ the in-going wave is exponentially dumped by the same reason. If we require 
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that the solution is decreasing in sector 3 of the Riemannian surface Sp and therefore the 
in-going and outgoing waves have equal amplitudes on the line -R+, we conclude that the 
in-going wave in T_ region and out-going wave in T+ region both have nonzero amplitudes. 
In quasiclassical expression for the wave function we could not retain the exponentially small 
part of solution simultaneously with exponentially large and therefore could not notice the 
outgoing wave in r_ region and the in-going wave in T+ region. 

Let us now analyze the requirement ([76iB). The solution of the equation (^5|) or (|^) is 
dehned on the Riemannian surface Sp (which is topologically two-dimensional sphere) with 
six punctured points: pi, p 2 , 0,1 - the turning points and two inhnities in S± components of 
Sp which are singular points. There exist six nontrivial basic cycles 71 ,..., yg around each of 
the punctured points on Sp. The procedure of construction of the quasiclassical solutions 
in complex domain gives the wave functions with trivial monodromy along the basic cycles 
7 i, ..74 around all the turning points 0, l,pi,p 2 - So we should only require that the wave 
function transforms into itself when prolonged along the pass 75 or 75 around the inhnities 
in S± components of Sp. It suffices to satisfy the condition only on one of the cycles because 
the other one is 


76 = (71 • 72 • 73 • 74 • 75) ^ (84) 

(bullet denotes the composition law in fundamental group) It will be convenient to consider 
the cycle T shown on Fig. (one could easily check that it is not a composition of cycles 
7 i, ...,74 around the turning points). One could see that we could cover all the pass T by 
two canonical regions = 21J5U3 and ^4 = 41J5U3 (see Fig. |^). Let us prolong the 
solution from the region 3 containing the R- line and the point A G F (Fig. ^) corresponding 
to the value t = 0 and t = 1 of the parameter along F to the region 5 containing Vp line 
and the point B eT corresponding to f = 1/2. 

The canonical region contains two Stokes lines I 5 and Ig. Therefore, there exists two 
different fundamental systems of solutions in this region ^'‘12 and $ 1 * 2 - The solution which 
decrease in the sector 3 of Rj? is 
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In order to clear up the behavior of this solution in sector 5 which belongs to the same 
canonical region D/^ as the sector 3 we express the solution 'h through the fundamental 
system of solutions dehned by the Stokes line /§ in canonical region D 4 . The transition 
matrix is 


^ 3,3 — 


exp I / PsdS 


exp / PsdS 


Here the integral is taken along the interval (0,1) corresponding to T_ line on the Riemannian 
surface Sp- So the wave function has the form 


T = exp J PsdS 


on the line V+ where $2 is increasing when p ^ —00 along 

Let us consider the form of wave function in the region D 4 ,. There are two Stokes lines 
and Ig in this region. In terms of fundamental system of solutions dehned by the Stokes 
line I 5 the wave function has the form: 

T = 04*1*]® T 64 $2^, p G D4 


When expressed through the fundamental system of solutions dehned by the Stokes line Ig 
the wave function becomes 


^ = exp - J PsdS 


where has the property +cxd when p —> — cx) in V+ 


region. 



From the symmetry of the construction with respect to real line we have 




( 90 ) 


SO the two expressions for the wave function 


exp 


and 


agree if 



PJS = 1 


( 91 ) 


Taking in account the fact that the wave function is in-going wave on T_ line and outgoing on 
T+ line, as it is explained earlier, we see that the integral in the argument of the exponent 
could be written as the integral over the closed path which goes around the cut made 
from the point p = 0 to p = 1 on both of the Riemannian spheres S± constituting the 
Riemannian surface Sp- So the quasiclassical wave function satisfy the requirement ([76|,B) 
if the quantization condition 

^ = KkeZ (92) 


holds. 

Collecting the results obtained in this section we conclude that in zero order on h the 
quasiclassical solution of equation (p5D dehned on the Riemannian surface Sp exists if the 

) hold. Taking in account the explicit expression 


two quantization conditions (|80D and ( 

(^) for Pg and calculating the integrals entering the quantization conditions we obtain the 
equations dehning the quasiclassical spectrum for the bound motion of the selfgravitating 
dust shell: 


Pi 


/ n JO 1 M 3 , f2m\ 

J 2 8m^ 2 \ M) ]\/[2 

V m? 


2 




PsdS = - 






2 m 


1 

2 


f m\ 

71 + Sarccos — 

\Mj 


~l 2 






777 ^ I \ 

——arccos — 

Mj 


( 93 ) 


Denoting cos a = m/M we obtain 
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ke z 


( 94 ) 


i f 1 = i — -—^-In (2 cos a) + a ctg 2 a = TiCk 

^ 4 cos a 

= + =<(^ + 2) 

The behavior of the functions standing in the left hand side of this equation is shown on 
Fig. Large k and large positive n correspond to wormhole states while large negative n 
and finite k correspond to the black hole states as one could see. 

Let us consider the spectrum in the black hole case. This corresponds to the values 
m < M < 2m of parameters on Fig. In this case the function /i changes in hnite 
limits from /i = 7r/(3\/3) — 1/2 at m/M = 1/2 to fi = 3/4 — ln2 at m/M = 1. So if the 
Schwarzschild mass m is given, then k could change only in finite limits 




In 2 


(95) 


So the black hole state with given Schwarzschild mass m (the only parameter that an observer 
at inhnity could measure) is a superposition of the states with different k and is in fact a 
mixed state, having nonzero entropy. Besides, we have the inequality 


m > mpi/, 


(96) 


\\6V3 dvr/ 

for the quasiclassical black hole spectrum because otherwise the inequality (p5| ) have no solu¬ 
tions. This means that the minimal black hole mass exists (if we suppose that quasiclassical 
spectrum is valid for low energy values). 

The mass spectrum for the black hole model under consideration was found in [Q in the 


large black hole limit using another technic. Let us compare the two spectra. From (|94D we 
hnd 
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= —{3k -|- 2n -|- 1) 


( 97 ) 




We see that although each of quantization conditions found in does not hold for the 
quasiclassical spectrum, their combination (0 does hold. 
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Let us consider the limit M > m, m ^ M which corresponds to the case when the 
turning point of the classical motion of the shell pi —>■ oo. In this limit a —>■ 0 and the last 
item in the second of equations (^) becomes large compared to the others. Then the second 
quantization condition (H) takes the form in this limit 




m2 

= n + 


(98) 




which coincides qualitatively with the quantization condition of |Q. So the two spectra 
coincide in the limit m —> M. 

VI. THE QUASICLASSICAL SPECTRUM FOR THE COLLAPSING SHELL. 

In the previous section we obtained the discrete mass spectrum for the bound motion of 
the selfgravitating shell. One of the quantization conditions (^,A) appeared to be similar 
to the Bohr-Sommerfeld quantization condition on the trajectories of the bound motion of 
mechanical system in nonrelativistic quantum mechanics. The other quantization condition 
(0 ,B) has a different origin being the requirement for the wave function to be a regular 


function on the Riemannian surface Sp- The first requirement (pbiA) will disappear if 
we consider the situation of gravitational collapse when the classical trajectory goes from 
p = oo to p = 0 (see Fig.|I|c). The wave function is the superposition of in-going and out¬ 
going waves near oo and does not decrease in this region. So the hrst of the quantization 
conditions (|3), (0) disappear. The principally new feature of our model compared to the 
ordinary nonrelativistic quantum mechanics is that the second condition (|^ which stems 
from the requirement (|76| ,B) for the wave function to be an unambiguous function on the 
Riemannian surface Sp does not disappear in the case of unbound motion and gives the 
discrete spectrum for the mass of the collapsing shell. 

The collapse case corresponds to the values of parameters m > M. In this case the 
integral entering the quantization condition (|9^) takes the form 
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(99) 


Denoting ch(a) = m/M we obtain 



( 100 ) 


The graph of this function is presented on Fig. We see that when the Schwarzchild 
mass is fixed there are only finite number of quantum states of the collapsing shell, because 
fi changes from 3/4 — ln2 to 1/2 for all the values 1 < m/M < oo. So 



( 101 ) 


The spectrum is continuous as far as it depends on a continuous parameter M but each level 
m = const is finite degenerate similarly to the black hole case. Similarly to the black hole 
case we see that equation ( p.01|) has solutions only if 


m > ^/nmpl 


( 102 ) 


which means that the minimal mass exists for the states of the collapsing shell in the 
quasiclassical spectrum. 

Let us consider now the collapsing null-dust shell ( M = 0 ). We write according to the 
Dirac quantization procedure the quantum version of the constraint (0). It is the following 
finite difference equation: 



( 103 ) 


We look for the quasiclassical solutions on the whole Sp of this equation. The corresponding 
Hamilton-Jacobi equation 



( 104 ) 


has the solutions 
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F, = ±ln 


(106) 


The integral entering the quantization condition (^) is 


1 


/ 



(106) 


0 


and (P^) gives the spectrum 


= Tirripik, k E Z 


(107) 


proposed by Bekenstein, Mukhanov and other authors. 

Let us consider the quasiclassical wave function in the case m > M. We would expect 
that the quasiclassical wave function will be concentrated near the classical trajectory of 
the collapsing shell in the phase space. Indeed, the in-going wave quasiclassical solution 
corresponds to this trajectory. But it turns out that we could not set the amplitude of the 
out-going wave to be equal to zero. We require that the wave function in R_ region to be 
equal to the decreasing solution of (p^ (otherwise it grows inhnitely with p —>■ cxo in 
region). Then we prolong of the solution to the whole Riemannian surface Sp following 
the procedure described in sections El and 0. We hnd that the in- going and out-going 
waves enter the wave function with equal amplitudes on Rp line. Therefore we see that the 
situation resembles the reflection from a potential wall in nonrelativistic quantum mechanics. 
The horizon point p = 1 is a turning point = 0 for the quasiclassical motion if we choose 
the regular coordinate u on the Riemannian surface Sp near the horizon. So the equation 
(^) describes the situation when the stationary in-going wave reflects completely from the 
horizon. In T_ region the wave function is an in-going wave and in T+ it is an out-going 
wave similarly to the case of bound motion. 


VII. CONCLUSION. 


In this paper we considered the quantum mechanical model of the black hole that consists 
of the selfgravitating thin dust shell. The Schroedinger equation for this model is a hnite 
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difference equation (|^) with the finite shift of the argument of the wave function along the 
imaginary axis |^. Therefore, the equation must be considered on a Riemannian surface Sp 
where the function (|^ is regular. The analysis of the equation on this Riemannian 
surface gives the mass spectrum for the black hole model in the quasiclassical approximation. 
While constructing the quasiclassical solutions on a complex manifold we need to determine 
the approximate solutions of the differential equation in different canonical domains of the 
complex manifold. Then we glue the global solution from the solutions dehned in different 
canonical regions with the help of transitions matrices between these regions. In the section 
0this method was used in order to hnd the quasiclassical solutions of the equation (^). The 
quasiclassical solution satisfying the two requirements (|^,A,B) exists if the coefficients of the 
equation satisfy two quantization conditions (|80D and (^). One of the obtained quantization 
conditions (^) follows from the boundary conditions at infinities in R± regions of Kruskal 
space-time ( requirement ([76|,A)) and is obtained by the same procedure as Bohr-Sommerfeld 
condition of nonrelativistic quantum mechanics. But the other (^) which follows from 
(0 ,B) and is expressed as the requirement of the trivial monodromy of the wave function 
along nontrivial cycle T (Fig. ^ on the Riemannian surface Sp with punctured singular 
points of the equation (P^l). It has no analogy in ordinary quantum mechanics. 

Being combined together the quantization conditions (|94D dehne the discrete mass spec¬ 
trum of the model, which depends on two quantum numbers. Furthermore, one of the 
quantization conditions (|9^) remains valid for the unbound motion where it takes the form 
( |100|) . The spectrum of unbound motion is continuous but each level is hnite degenerate. 
The black hole spectrum of both bound and unbound motions are bounded from below 
m > ^/nmpl. For the null-dust M = 0 the remaining quantization condition gives the dis¬ 
crete spectrum of unbound motion which turns out to be the spectrum found by Bekenstein 
and Mukhanov IHI. 
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FIGURES 

FIG. 1. The space-time with selfgravitating thin shell. The black hole case (a) the wormhole 


case (b) and the case of collapsing shell (c). 

FIG. 2. The structure of Stokes lines for the potential with two turning points zq and zi. 

FIG. 3. The real section of the Riemannian surface Sp is a cross naturally identified with the 
configuration space of classical dynamical system. 

FIG. 4. The Stokes lines on Sp near the turning points pi and p 2 - 

FIG. 5. The Stokes lines near the singular points p = oo in S±. 

FIG. 6. The Stokes lines in the vicinity of the horizon p = 1 and near the singularity p = 0. 

FIG. 7. The global structure of Stokes lines on the Riemannian surface Sp. 

FIG. 8. Functions defining the mass spectrum. 
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